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PREFACE 

This  topical  report  was  prepared  by  Science  Applications,  Inc. 
(SAI)  in  compliance  with  the  requirements  of  Defense  Nuclear  Agency 
Contract  DNA001-76-C-0023.  The  program  was  sponsored  by  the  Strategic 
Structures  Division  of  the  Shock  Physics  Directorate  of  DNA  and  was 
performed  during  the  period  01  October  1975  through  30  September  1976. 
The  technical  representative  of  the  sponsor  was  Capt.  Jerry  R.  Stockton. 

The  program  at  SAI  was  performed  under  the  direction  of 
Dr.  William  R.  Seebaugh.  He  was  assisted  in  the  computer  analyses  by 
Mr.  Elwood  E.  Zimmerman  and  Ms.  Susan  J.  Rose. 
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Section  1 


INTRODUCTION 

The  crater  ejecta  environment  for  high  explosive  (HE)  and 
nuclear  bursts  has  been  under  investigation  by  Science  Applications, 

Inc.  (SAI)  for  several  years.  Earlier  work  on  the  dynamic  crater 
ejecta  model,  supported  by  the  Defense  Nuclear  Agency  (DNA),  was 
conducted  under  Contracts  DNA001-73-C-0088,  DNA001-74-C-0102  and 
DNA001-75-C-0102  during  the  period  20  November  1972  through  30  October 
1975.  The  basic  ejecta  model  was  developed  during  this  period.  The 
present  effort,  which  was  performed  under  Contract  DNA001-76-C-0023, 
was  concerned  with  updating  the  model  to  incorporate  the  results  of 
the  most  recent  theoretical  calculations  and  experimental  investigations. 
The  period  of  performance  of  this  contract  was  01  October  1975  through 
30  September  1976.  This  report  describes  the  current  ejecta  model 
and  thus  provides  a summary  of  the  entire  four-year  effort. 

The  ejecta  mass  represents  a substantial  fraction  of  the  total 
mass  excavated  during  the  formation  of  an  explosive  crater.  This  mass 
has  been  determined  for  experimental  craters  by  measuring  the  thicknesses 
of  the  post-event  ejecta  blankets,  A number  of  empirical  models  for 
the  thickness  of  the  ejecta  blanket  have  been  developed;  see,  for 
example.  Post  (Reference  1)  and  McGetchin  et  al . (Reference  2).  These 
models  are  based  upon  measurements  obtained  primarily  from  a large 
number  of  HE  and  nuclear  craters.  Limited  data  are  also  available 
for  the  Arizona  Meteor  Crater  and  estimates  of  the  ejecta  thickness  at 
the  crater  rim  have  been  made  for  several  lunar  craters. 

An  alternative  method  of  obtaining  estimates  of  the  ejecta 
blanket  thickness  considers  the  entire  history  of  the  crater  formation 
and  ejecta  deposition.  In  this  approach,  which  is  the  subject  of  this 
report,  the  crater  formation  and  ejecta  throwout  processes  are  modeled, 
the  trajectories  of  representative  ejecta  fragments  are  calculated 
and  the  ejecta  blanket  thickness  is  obtained  by  summing  the  contributions 
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of  the  fragments  impacting  the  surface.  The  procedure  also  provides 
a description  of  the  airborne  ejecta  cloud  prior  to  deposition  of 
the  ejecta  fragments  on  the  surface.  The  problem  is  delineated 
schematically  in  Figure  1.  The  left  side  of  the  sketch  illustrates 
the  positions  of  the  representative  fragments  a few  seconds  after  the 
burst.  These  fragments  impact  the  surface  within  about  a minute  of 
the  burst,  forming  the  ejecta  blanket  as  shown  on  the  right  side  of 
the  sketch.  The  inner  region  of  the  ejecta  blanket  is  characterized 
by  a continuous  distribution  of  material  that  completely  covers  the 
original  ground  surface.  At  larger  ranges,  there  is  insufficient 
ejecta  mass  to  cover  the  surface  and  the  distribution  is  discontinuous. 
The  ejecta  mass  is  defined  as  the  amount  of  crater  material  that 
crosses  the  original  ground  surface  and  impacts  beyond  the  apparent 
crater  radius  (R.  in  Figure  1);  this  definition  excludes  both  fallback 

a 

and  upthrust.  The  crater  lip  radius  is  denoted  by  Rg-|. 


Section  2 


GENERAL  DESCRIPTION  OF  EJECTA  MODEL 

The  ejecta  model  is  based  upon  recent  crater  calculations, 
selected  low-yield  ejecta  data  (primarily  from  HE  events)  and  current 
models  of  the  buoyant  fireball  and  the  associated  wind  field.  The 
model  is  illustrated  in  flow  chart  form  in  Figure  2.  The  inputs  to 
the  model  are: 

(1)  Theoretical  velocity  distribution  for 
continuous  hydrodynamic  material  emerging 
from  the  crater 

(2)  Empirical  models  for  ejecta  mass,  fragment 
size  distribution  and  maximum  fragment  size 

(3)  Experimental  observations  of  fragment  velocity 
for  high  explosive  (HE)  shots 

(4)  Vortex  flow  model  (VORDUM®)  for  the  fireball 
and  winds  which  affect  ejecta  fragment 
trajectories 

(5)  Compressible  flow  aerodynamic  drag  model. 

Items  (1)  through  (3)  above  define  the  crater  source.  Ejecta  fragments 
produced  by  this  source  are  transported  through  the  flowfield  defined 
by  the  VORDUM  program  (Item  (4)  above]  using  the  aerodynamic  drag 
model  for  ejected  crater  material  (Item  (5)  aoove].  The  source  model 
and  the  trajectory  calculation  form  the  complete  ejecta  model  as^ shown 
in  Figure  2.  The  method  of  analysis  is  applicable  to  the  modeling  of 
both  the  ground  deposition  of  crater  ejecta  and  the  airborne  ejecta 
cloud.  The  deposition  results  are  obtained  in  the  form  of  the  cumulative 
mass  and  number  of  fragments  impacting  the  ground  plane  per  unit  area, 
the  ejecta  depth,  the  fragment  size  class  (minimum  and  maximum  fragment 
diameters)  and  the  minimum  and  maximum  values  of  the  impact  time, 

®V0RDUM  is  an  acronym  for  a vortex  ^st  model  which  was  developed  to  describe 
tho  fluid  and  particulate  motion  associated  with  the  rise  of  a buoyant 
fireball  through  the  earth's  atmosphere. 
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velocity,  angle,  momentum  and  kinetic  energy  as  functions  of  range 
from  the  burst  point.  The  ejecta  cloud  environment  is  produced 
simultaneously.  During  the  VORDUM  trajectory  calculations  (Figure  2), 
trajectory  parameters  in  the  form  of  range  and  altitude  versus  time 
are  recorded  for  each  representative  fragment.  This  information  is 
combined  with  the  crater  source  description  to  give  airborne  ejecta 
fragment  mass  and  number  densities. 


Section  3 


THE  EJECTA  SOURCE 

The  most  impc  ^'tant  assumptions  of  the  ejecta  model  are  related 
to  the  ejecta  source.  The  source  description  must  include  the  mass 
ejection  rate,  the  ejection  velocity  (speed  and  angle)  and  the  fragment 
size  distribution  as  functions  of  radius  and  time.  Empirical  correlations 
have  been  derived  for  the  total  ejecta  mass  (Reference  3)  and  for  the 
relative  origins  of  ejecta  fragments  within  the  crater  (Davis  and  Carnes, 
Reference  4).  The  latter  correlation  is  shown  in  Figure  3 for  event 
MINERAL  ROCK,  a 100-ton  TNT  shot  on  granite.  As  shown,  the  ejecta 
missiles  originate  from  the  upper  portion  of  the  crater  away  from  the 
center  of  the  charge.  Very  little  ejecta  comes  from  the  region  near 
the  axis.  Although  cratering  studies  have  produced  correlations  such 
as  those  shown  in  Figure  3,  attempts  to  construct  an  ejecta  source 
model  completely  from  experimental  data  have  been  unsuccessful.  Data 
do  not  exist  on  the  rate  of  ejection  of  material  from  the  crater.  The 
in-flight  photography  of  ejecta  fragments  may  be  analyzed  to  give 
representative  fragment  ejection  velocities,  but  such  analyses  do  not 
relate  mass  to  velocity.  It  is,  therefore,  necessary  to  incorporate 
results  of  theoretical  calculations  into  the  ejecta  source  model  to 
complete  the  source  description. 

A number  of  two-dimensional  hydrodynamic  calculations  using 
advanced  elastic-plastic  material  models  have  been  performed  recently. 

The  methods  have  not  been  generally  successful  in  predicting  the  final 
crater  dimensions.  It  appears  that  the  poor  agreement  for  the  final 
crater  volume  is  a result  of  failure  to  properly  predict  the  late- 
time motion,  which  occurs  at  low  velocities.  The  mass  of  ejecta 
leaving  the  crater  with  higher  velocities  is  not  as  dependent  on  the 
material  model  as  the  low-velocity  ejecta,  and  consistent  results  have 
been  obtained  for  the  high-velocity  material  by  assuming  that  the  ejecta 
motion  is  driven  by  the  amount  of  energy  coupled  to  the  ground  (defined 
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Figure  3.  Zones  of  Origin  of  Crater  Ejecta 


as  the  total  energy  below  the  ground  surface  at  1 ysec  after  burst). 

The  results  of  an  analysis  of  the  properties  of  the  crater  material 
crossing  the  original  ground  surface  for  a number  of  calculations  are 
presented  in  Figure  4 in  the  form  of  average  vertical  velocity  (defined 
as  the  ratio  of  the  momentum  flux  to  the  mass  flux)  ver-sus  the  cumulative 
ejecta  mass  with  greater  average  vertical  velocity.  Curve  1 represents 
the  average  relationship  for  high  strength  materials  (wet  tuff  and 
granite),  whereas  Curve  2 is  the  average  for  two  calculations  of  the 
nuclear  event  JOHNIE  BOY  (0.5  kt  yield  at  a depth  of  burst  of  0.6  m in 
alluvium).  These  curves  suggest  the  source  relationship  given  by  the 
dashed  line  in  Figure  4 as  an  upper  velocity  bound  for  a nuclear 
surface  burst.  The  equation  of  this  line  may  be  represented  by  the 
function 


aM 

W 


^min.z  y-2 
sine 


(1) 


where  3M/3V  is  the  ratio  of  the  differential  ejecta  mass  6M  associated 

with  differential  ejection  velocity  6V,  is  the  ejecta  mass,  6 is  the 

ejection  angle  relative  to  the  horizontal  and  V is  the  ejection  velocity. 

This  function  is  normalized  to  give  the  empirically  determined  ejecta 

mass  (see  below)  at  a minimum  value  of  the  vertical  component  of  the 

ejection  velocity  V , of  6.4  m/sec.  The  minimum  velocity  considered 
min,z 

is  then 

V _ (^min,  z)  _ 6.4  m/sec  /p\ 

^min  " sin0  sine  ' ' 


The  total  ejecta  mass  is  determined  from  correlations  of 
experimental  cratering  data.  Layson  (Reference  3)  has  quoted  the 
relationship 

= 2.4  (3) 

for  the  apparent  crater  mass  (M, , megatons)  for  megaton-range  surface 

d 

bursts  on  rock,  where  W is  the  explosive  yield  in  megatons.  Experimental 
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Figure  4.  Cumulative  Ejecta  Mass  From  Crater  Calculations 


data  (Reference  3)  also  suggest  that  the  ejecta  mass  (in  megatons) 
is  given  by 


= 0.6  M,  (4) 

e a 

or,  combining  Equations  (3)  and  (4), 

= 1.44  . (5) 

e 

The  ejecta  mass  includes  all  excavated  material  crossing  the  original 
ground  surface,  a large  part  of  which  emerges  from  the  crater  at  rela- 
tively low  velocities  and  comes  to  rest  on  the  crater  lip  (Figure  1). 
Material  with  initial  vertical  velocity  below  6.4  m/sec  is  not  included 
in  the  ejecta  mass  M^. 

The  apparent  crater  mass  given  by  Equation  (3)  corresponds  to 

3 

a cratering  efficiency  of  0.84  m /ton;  this  value  is  in  agreement  with 
the  mean  cratering  efficiency  given  by  Post  (Reference  1)  for  hard  rock. 
Liner  (Reference  5)  examined  the  geology  and  height-of-burst  dependence 
of  the  cratering  efficiency  and  derived  the  following  relationship  for 
apparent  crater  mass  as  a function  of  yield  and  scaled  height  of  burst 
(SHOB): 


M^  = a^  p W KgK^  exp(-d2  SHOB)  (6) 

where  p is  the  in  situ  density  of  the  cratering  medium  (kg/m  ),  Kg  and 
Kg  are  geology  and  source  correction  terms  and  a-|  and  a2  are  given  by 


Source 

Nuclear  (LYTW 
and  HYTW) 

HE 


4.8  X 10"^ 
2.06  X 10" 

4.8  X 10'^ 
2.06  X 10' 


SHOB  Range  (m/kt^^^) 


1.007  -1.56  to  1.1 

0.069  -19.5  to  -1.56 

0.325  -1.79  to  5.18 

0.069  -19.5  to  -1.79 


The  cratering  medium  density  and  the  geology  correction  term  were 
determined  by  classifying  HE  data  by  generic  geologies: 
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Generic  Geology 

p(kg/m^) 

Kg 

Hard  Rock 

2600 

0.5  _+  0.2 

Dry  Sand-Sandstone 

1500 

0.58  + 0.2 

Alluvium 

1600 

1.00  + 0.23 

Clay-Shale 

1760 

1.89  + 0.89 

Saturated  Soil 

2100 

4.41  + 0.64 

The  source  correction  term  distinguishes  HE,  low  yield-to-weight  (LYTW) 
nuclear  and  high  yield-to-weight  (HYTW)  nuclear  sources: 


HE: 

^s  = 1 

^ 5.71 

(-1.79  m/kt^''^<SH0B<5.18  m/kt^''^) 

1 

[2.1 

(-19.5  m/kt^''^<SH0B<-1.79  m/kt^'^^; 

LYTW 

Nuclear: 

's  = 

1 

HYTW 

rO.25 

(-0.3  m/kt^''^<SH0B<l.l  m/kt^^^ 

Nuclear: 

S = 

1 0.25  - 

0.18(1+3.3  SHOB)  (-1.56  m/kt^''^  < SHOB  <-0.3  m/kt^''^ 

1 

ll 

(SHOB  <-1.56  m/kt^''^) 

Equation  (4)  is  applied  to  determine  the  ejecta  mass. 

Equation  (6)  gives  a cratering  efficiency  of  0.55  m^/ton  for 
hard  rock  at  SHOB  = 0.  This  is  about  35  percent  below  the  aforementioned 
value  quoted  by  Layson  (Reference  3)  and  Post  (Reference  1)  and  is  at  the 
lower  end  of  the  scatter  band  given  by  Post.  Both  Liner  (Reference  5) 
and  Post  (Reference  1)  scaled  the  apparent  crater  mass  linearly  with 
yield,  whereas  Layson  (Reference  3)  scaled  mass  with  yield  to  the  0.9 
power.  The  effects  of  these  differences  are  within  the  scatter  of  the 
data  base  from  which  the  various  correlations  were  derived.  The 


approach  developed  by  Liner  [Equation  (6)]  has  been  incorporated  into 
the  current  ejecta  model.  For  the  case  discussed  in  Section  5,  a 
cratering  efficiency  of  0.84  m /ton  [Equation  (3)]  was  used  for  a 
1 MT  surface  burst  over  hard  rock. 


The  distribution  of  the  crater  mass  over  the  range  of  possible 
ejecta  fragment  sizes  is  an  important  ingredient  of  the  ejecta  source 
model.  Studies  of  cratering  in  rock  media  such  as  basalt  and  granite 
(Reference  3)  indicate  that  the  size  distribution  has  the  form 


3d 


0.5  M 
e 

m 


(7) 


where  8M/3a  is  the  ratio  of  the  differential  mass  6M  associated  with 

the  differential  fragment  diameter  6a,  over  a range  of  fragment  diameter 
-4 

from  about  10  to  1 m.  The  term  a^^^  is  the  diameter  of  the  largest 
ejecta  fragment.  Gault,  et  al.  (Reference  6)  compiled  data  for 
explosive  and  impact  cratering  events  and  observed  that  the  correlation 
for  the  mass  of  the  largest  fragment 

m^^  = 6.6  X 10^  (8) 


was  a reasonable  representation  of  the  data  for  a variation  of  14  orders 
of  magnitude  of  M^,  The  constant  of  proportionality  gives  the  maximum 
fragment  mass  m^  in  kg  when  is  expressed  in  megatons.  Assuming 
spherical  fragments  the  maximum  fragment  diameter  a^(»n)  is  give  by 


The  in  situ  density  is  used  in  the  current  model.  The  size  distribution 
given  by  Equation  (7),  with  the  upper  limit  given  by  a^j^,  is  assumed 
to  apply  to  all  cratering  media  (hard  rock,  sandstone  and  so  forth) 
except  those  clearly  defined  with  alluvial  basins  or  sands  with  water 
table  levels  below  the  apparent  crater  depth. 
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The  size  distribution  for  incohesive  soils  (alluvium  and  dry 
sand)  is  more  accurately  described  by  the  relationship  derived  for 
Nevada  Test  Site  alluvium  (Reference  3) 

= (10) 

9a  a 

The  maximum  ejecta  fragment  diameter  for  such  incohesive  soils  is  likely 
to  be  more  closely  associated  with  the  maximum  in  situ  diameter  than 
with  the  ejecta  mass  [Equation  '8)].  For  these  situations,  the  in  situ 
maximum  diameter  is  employed  as  the  upper  limit  to  the  distribution 
given  by  Equation  (10). 

There  is  currently  no  experimentally  based  method  for  selecting 
size  distributions  between  the  limits  represented  by  Equations  (8)  and 
(10).  Fortunately,  the  ejecta  environment  is  not  extremely  sensitive 
to  this  parameter.  The  effects  of  using  the  in  situ  maximum  diameter 
in  place  of  the  correlation  of  Equation  (8)  (for  alluvium  and  dry  sand] 
increase  as  the  difference  between  the  respective  maximum  diameters 
increases.  This  is  a result  of  the  influence  of  air  drag.  Further 
research  is  required  to  more  accurately  define  these  effects  for 
incohesive  cratering  media. 

The  crater  mass-velocity  relationships  derived  from  the  cratering 
calculations  and  the  experimentally  derived  fragment  size  distributi'>ns 
cannot  be  truly  coupled  because  the  process  of  comminution  of  the 
crater  material  has  not  been  described  either  theoretically  or 
experimentally.  The  approach  taken  in  the  current  ejecta  model  is  to 
distribute  the  ejecta  mass  over  the  velocity-size  spectrum  using  the 
aforementioned  relationships,  and  to  exclude  combinations  of  velocity 
and  size  that  appear  unreasonable  based  upon  an  analysis  of  additional 
experimental  data. 


The  basis  for  a limiting  velocity-size  relationship  is  provided 
by  the  results  of  an  extensive  in-flight  ejecta  analysis  using  photography 
of  HE  events  (Reference  7).  The  analysis  gives  the  sizes  and  ejection 
conditions  for  a large  number  of  large  ejecta  fragments.  The  results 
for  event  MIDDLE  GUST  III,  a 100-ton  surface  tangent  TNT  shot  on  a 
layered  (clay  over  shale)  medium,  are  presented  in  Figure  5®.  The 
fragment  diameters  were  obtained  by  assuming  spherical  fragments.  The 
data  are  bounded  by  two  limits: 

"max  = 566 


and 


V = C V 
m max  a 


m 


(12) 


where  velocities  are  in  m/sec.  The  material  ejected  at  the  highest 
hydrodynamic  velocities  near  the  burst  (Figure  4)  must  either  be 
vaporized  or  pulverized  by  the  shocks  measured  in  hundreds  of  kilobars. 
It  is  not  reasonable  to  associate  these  high  velocities  with  fragments 
of  significant  size.  The  absolute  maximum  velocity.  Equation  (11), 
performs  the  required  limiting  function,  with  the  value  of  the  constant 
determined  from  the  data  of  Figure  5.  The  1/6-power  dependence  of 
velocity  on  yield  was  suggested  by  Sauer  and  Cooper  (Reference  8)  and 
observed  experimentally  by  Wisotski  (Reference  7).  Qualitative 
observations  of  results  of  field  surveys  and  aerial  photography  for  HE 
and  nuclear  events  indicate  that  the  largest  ejecta  fragments  do  not 
reach  the  largest  ranges  from  the  burst.  The  use  of  a size-dependent 
maximum  velocity.  Equation  (12),  produces  the  observed  results.  The 
form  of  this  limit  was  also  obtained  from  the  data  of  Figure  5. 


small  number  of  fragments  (9  out  of  a total  of  151)  which  required 
large  corrections  to  obtain  the  true  azimuthal  planes  of  the  trajectories 
were  eliminated  from  the  sample.  See  Reference  7 for  an  explanation 
of  the  correction  procedure. 
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Figure  5.  Correlation  of  Ejecta  Velocity  and  Fragment  Diameter 


The  parameter  C is  a function  of  the  explosive  type  (HE  or  nuclear) 
and  the  depth  of  burst.  Tentative  values  for  this  parameter  are  0.5 
for  surface  tangent  HE  (for  example,  event  MIDDLE  GUST  III),  0.3  for 
half-buried  HE  and  0.1  for  megaton-range  nuclear  surface  bursts. 

Equation  (12)  is  required  only  for  cunesive  media;  the  limit  provided 
by  Equation  (11)  is  sufficie^c  to  bound  the  fragment  velocities  for 
incohesive  media  for  which  the  maximum  fragment  diameter  equals  the 
in  situ  value.  Note  that  the  size-dependent  velocity  limit  shown  in 
Figure  5 was  calculated  using  the  experimentally  determined  ejecta 
mass  in  Equation  (8). 

The  relationships  given  in  Equations  (1)  through  (12)  define 
the  overall  characteristics  of  the  ejecta  source.  In  order  to  complete 
the  source  model,  the  points  of  origin  of  the  ejecta  fragments  within 
the  crater  must  be  specified.  The  experimental  ejecta  origin  description 
illustrated  in  Figure  3 is  modeled  effectively  by  a crater  model 
developed  by  Rosenblatt  and  Eggum  (Reference  9).  This  model  simulates 
the  time  phased  development  of  the  crater  as  described  by  the  same 
hydrodynamic  cratering  calculations  used  to  determine  the  mass -velocity 
relationships  shown  in  Figure  4.  The  ejection  velocity  and  mass  flux 
across  the  original  ground  surface  are  specified  as  functions  of  time 
(consistent  with  the  dashed  line  in  Figure  4).  For  cohesive  media, 
each  increment  of  mass  is  partitioned  into  a range  of  fragment  sizes 
according  to  Equation  (7),  consistent  with  the  limits  specified  by 
Equations  (11)  and  (12).  Equations  (10)  and  (11)  perform  the  same 
functions  for  incohesive  soils.  The  result  is  a matrix  of  representative 
ejecta  fragments,  with  each  fragment  characterized  by  its  ejection 
velocity  (speed  and  angle),  diameter,  ejection  time. location  within  the 
developing  crater  and  the  ejecta  mass  that  it  represents.  This  matrix 
provides  the  initial  conditions  for  the  ejecta  trajectory  calculations. 


Section  4 


EJECTA  TRAJECTORIES 

The  primary  effects  of  the  explosion  are  tc  form  the  crater,  to 
generate  a hot  low  density  region  near  the  burst  point  and  to  create 
strong  afterwinds  which  perturb  the  ejecta  trajectories.  It  is  of 
interest  to  know  more  of  the  motion  through  the  very  early  fireball 
environment,  but  it  is  beyond  the  -:tate  of  the  art  to  describe  the 
details  of  the  interaction  of  the  ejected  earth  material  and  the 
developing  early  fireball.  The  most  reliable  way  to  circumvent  this 
diff-'culty  is  to  make  good  calculations  of  the  trajectories  after 
pressure  equilibration  while  adjusting  the  source  description  to  obtain 
good  agreement  with  the  available  HE  and  nuclear  test  data.  This  is 
the  approach  that  has  been  followed  in  the  development  of  the  current 
ejecta  model . 

A relatively  simple  hydrodynamic  model,  VORDUM,  has  been 
developed  to  describe  the  air  and  dust/debris  motion  inside  and  around 
a rising  cloud.  In  this  model,  the  early-time  dust  cloud  is  a buoyant 
vortex  whose  flowfields  entrain  dust  particles  and  distribute  them 
through  space  according  to  their  initial  position,  mass,  and  aerodynamic 
drag  characteristics.  Within  the  fireball  radius,  the  aerodynamic 
drag  is  calculated  for  air  at  one-tenth  ambient  density.  Trajectories 
of  larger  ejecta  fragments  are  perturbed  by  the  wind  field.  The 
boundary  condition  at  the  ground  — no  normal  velocity  component  at  the 
ground  surface  — is  satisfied  by  imposing  an  image  vortex  system  that 
moves  away  from  the  ground  in  the  negative  field.  As  the  cloud  rises, 
it  expands  adiabatically  and  also  mixes  with  the  ambient  air.  The 
rates  at  which  the  cloud  rises  and  expands  viere  obtained  from 
test  data  and  hydrodynamic  cloud  calculations. 
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A number  of  trace  fragment  trajectories  is  computed  for  each 
ejecta  environment  description  desired.  Cylindrical  syimietry  is  assumed 
at  all  times  during  the  calculation.  The  VORDDM  output  consists  of 
range  and  altitude  coordinates  describing  each  trace  fragment  trajectory. 

Each  trace  fragment  corresponds  to  many  actual  fragments;  the  ejecta 
mass  corresponding  to  each  representative  fragment  is  obtained  from 
Equations  (1)  and  (7)  for  cohesive  cratering  media,  and  from  Equations  (1) 
and  (10)  for  incohesive  soils. 

Perhaps  the  most  important  assumption  in  the  trajectory  calculation 
relates  to  the  drag  coefficient.  Previous  studies  have  used  values 
measured  for  smooth  spheres,  but  as  explained  below,  this  assumption  is 
not  appropriate  for  ejecta  fragments  and  the  present  model  employs  a 
different  drag  model.  The  standard  drag  curves  for  flow  over  a smooth 
sphere  are  given  in  Hoerner  (Reference  10).  The  range  of  Reynolds 
numbers^  of  interest  in  the  ejecta  trajectory  calculations  is  from  about 
10  to  10  . For  incompressible  flow,  the  drag  coefficient  remains 
nominally  constant  in  the  region  of  laminar  boundary- layer  flow  from 

*3  C 

Re  - 10  to  about  3 x 10  . In  this  regime,  a substantial  portion  of 
the  total  drag  is  the  pressure  drag  engendered  by  the  separation  of  the 
flow  away  from  the  sphere  near  its  equator.  At  a Reynolds  number  of 
about  3 X 10^  for  a smooth  sphere  the  boundary  layer  undergoes  transition 
to  a turbulent  state  and  the  separation  point  moves  aft  over  the  down- 
stream surface  of  the  sphere,  substantially  reducing  the  pressure  drag 
and  hence  the  total  drag.  Values  of  Cq  from  0.1  to  0.2  are  obtained 
at  Reynolds  numbers  greater  than  about  5 x 10^.  Compressibility  effects 
on  the  drag  coefficient  are  large  and  nonlinear  throughout  the  Reynolds 
number  range  of  interest.  The  influence  of  transition  to  a turbulent  boundary 


^Reynolds  number  is  the  ratio  of  inertial  to  viscous  forces:  Re  = pVa/ji, 
where  p = atmospheric  density,  V = fragment  velocity  relative  to  the 
atmosphere,  a = fragment  diameter  and  pi  = dynamic  viscosity.  Drag 
coefficient  is  the  ratio  of  the  drag  force  to  the  product  of  the 
dynamic  pressure  and  the  fragment  cross-sectional  area. 


Cn  = D/(l/2pv2)(ira2/4). 
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layer  is  evident  only  for  Mach  numbers^  below  the  critical  value  of 
about  O.f.  At  higher  Mach  numbers,  the  shock  waves  formed  in  the 
flow  near  the  equator  separate  the  boundary  layer  at  this  point  for 
all  Reynolds  numbers  of  interest  and  no  drag  reduction  occurs.  For 
Mach  numbers  greater  than  1.5,  the  drag  coefficient  is  approximately 
1.0,  independent  of  Reynolds  number. 

The  aforementioned  discussion  c'oplied  to  smooth  spheres  and 
spheres  roughened  to  about  3 percent  of  the  diameter  (about  the  roughness 
of  a golf  ball).  The  ejecta  fragments  under  consideration,  in 
particular,  the  larger  fragments,  are  not  spherical  and  the  values  for 
smooth  or  slightly  roughened  spheres  do  not  apply.  Drag  coefficients 
for  bluff  shapes  such  as  short  blunt  cylinders,  wedges  and  so  forth, 
with  lengths  nearly  equal  to  the  diameter  or  width,  range  from  about 

3 

0.8  to  1.2  for  Reynolds  numbers  greater  than  10  . Of  even  greater 
significance  than  the  generally  higher  level  for  bluff  shapes,  as  compared 
to  smooth  spheres,  is  the  observation  that  no  drag  decrease  occurs  in 
the  critical  Reynolds  number  range  for  bodies  with  sharp  edges  or 
blunt  aft  faces.  This  is  due  to  the  "fixing"  of  the  point  of  boundary- 
layer  separation  at  the  sharp  corner  of  the  body.  These  observations 
apply  to  data  obtained  for  cubes  tumbling  in  an  uncontrolled  manner 
(data  summarized  by  Hoerner,  Reference  10).  A value  of  Cq  of  0.75 
was  obtained  in  the  incompressible  range  for  cubes.  The  influence  of 
compressibility  on  the  drag  coefficients  for  bluff  shapes  and  cubes  is 
similar  to  that  for  smooth  spheres,  with  Cp  increasing  to  1.2  for  a 
Mach  number  of  about  1.4,  and  remaining  nearly  constant  at  that  value 
for  higher  Mach  numbers. 

Since  ejecta  fragments  have  been  observed  to  be  very  rough  and 
nonspherical  in  shape,  the  smooth  sphere  data  were  not  incorporated 
into  the  ejecta  model.  The  most  important  effect  is  the  elimination  of 
the  drag  decrease  at  the  critical  Reynolds  number.  The  simplest 
approach  employs  an  incompressible  arag  law  of  the  form 

® Mach  number  is  the  ratio  of  the  fragment  velocity  relative  to  the 
atmosphere  to  the  speed  cf  sound. 
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The  next  step  is  to  consider  the  effects  of  the  air  compressibility.  A 
cross-plot  of  the  data  given  by  Hoerner  (Reference  10),  delineating  the 
dependence  of  the  smooth  sphere  drag  coefficient  on  Mach  number  for 
Reynolds  numbers  of  interest  in  the  » jecta  study,  is  shown  in  Figure  6. 

o c 

The  range  of  drag  coefficients  giver  by  Equation  (13)  for  10  < Re < 3 x 10 
is  from  0.60  to  0.64;  thus,  the  incompressible  drag  coefficients  are 
too  low  even  for  smooth  spheres  above  a Mach  number  of  about  0.7. 
Calculations  are  currently  performed  using  a compressible  drag  law 
shown  by  the  dashed  lines  in  Figure  6.  This  compressible  drag  law  has 
the  following  characteristics: 

t It  obeys  Equation  (13)  for  incompressible 
flow,  defined  here  as  M<0.4  for  all  Reynolds 
numbers 

• It  fits  data  from  Hoerner  (Reference  10)  for 
supersonic  flow  over  bluff  objects  at  high 
Reynolds  numbers,  defined  here  as  M>1.2  for 
Re  >60 

# It  obeys  Equation  (13)  for  slow  viscous  flow 
at  all  Mach  numbers,  defined  here  as  Re  160. 

For  the  intermediate  ranges  of  Reynolds  and  Mach  numbers,  the  drag 
coefficient  is  determined  by  interpolating  linearly  between  Equation  (13) 
at  M = 0.4  and  = 1.2  at  M = 1.2.  This  procedure  is  shown  in  Figure  6 
for  two  Reynolds  numbers. 
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Figure  6.  Fragment  Drag  Coefficients 
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Section  5 


RESULTS  AND  DISCUSSION 

Calculations  of  the  ejecta  deposition  and  cloud  environments 
were  performed  for  a 1 MT  surface  burst  on  hard  rock.  A matrix  of 
12,000  trace  fragments  was  developed  using  Equations(l)  and  (7).  The 
fragment  diameter  interval  of  0.002  to  9 m was  divided  into  66  size 
classes  of  equal  mass.  The  absolute  maximum  velocity  was  566  m/sec 
[Equation  (ll)].  The  size-dependent  maximum  velocity  was  obtained 
from  Equation  (12)  using  C=0.1. 

The  predicted  ejecta  depth  is  shown  in  Figure  7 (current  model). 

The  distribution  begins  at  the  apparent  crater  radius,  0.14  km.  The 

peak  depth  of  3.2  m occurs  at  the  crater  lip  radius,  0.17  km.  The 

depth  then  decreases  with  increasing  range  out  to  the  maximum  range 

-4 

of  6.5  km.  The  initial  decrease  follows  approximately  an  R law.  The 

-2  5 

depth  between  about  1 and  3 km  range  decreases  by  about  R ' . The 

-5  -3 

distribution  at  larger  ranges  varies  as  R . The  R law  proposed  by 
McGetrhin  et  al . (Reference  2)  for  lunar  craters  is  also  shown  in 
Figure  7.  Finally,  results  obtained  using  the  exponential  model 
developed  by  Post  (Reference  1)  are  given  for  the  maximum  (95  percentile), 
median  (50  percentile)  and  minimum  (5  percentile)  ejecta  thicknesses 
using  a cratering  efficiency  of  0.84  m /ton  (equivalent  to  current  model). 
Post's  median  correlation  is  within  a factor  of  2 of  the  results  obtained 
using  the  current  model  out  to  a range  of  about  1.5  km.  Beyond  this 
range,  the  basic  characteristics  of  the  distributions  are  different 
(power  law  vs.  exponential)  and  Post's  median  correlation  drops  below 
the  results  obtained  using  the  current  model.  The  exponential  model 
thus  appears  to  severely  underpredict  the  ejecta  depth  at  large  ranges; 
however,  the  depth  in  this  discontinuous  ejecta  regime  is  not  a particularly 
significant  parameter  (see  Figure  8,  discussed  below). 
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The  maximum  (95  percentile)  and  minimum  (5  percentile)  curves 
of  the  Post  model  (Reference  1)  account  for  the  scatter  in  the 
experimental  data  base.  Tnis  scatter  is  primarily  a result  of  the 
azimuthal  variations  in  tfie  ejecta  distribution  pattern  commonly 
known  as  rays.  The  current  model  does  not  have  provisions  for  the 
determination  of  azimuthal  effects.  The  uncertainty  in  the  results 
introduced  by  model  uncertainties  other  than  azimuthal  variations  is 
estimated  to  be  about  a factor  of  3 relative  to  the  values  shown  in 
Figure  7.  The  estimation  of  95  and  5 percentile  limits  for  the  current 
model,  including  ray  effects,  merits  further  research  effort. 

The  ejecta  areal  mass  density  (mass  deposited  per  m surface 
area),  which  was  presented  in  previous  reports  as  the  principal  deposition 
parameter  (see  Reference  11,  for  example),  may  be  obtained  by  multiplying 
the  ejecta  depths  given  in  Figure  7 by  the  in  situ  density  of  the 
cratering  medium.  The  corresponding  areal  number  density  (number  of 
fragments  deposited  per  m surface  area)  is  presented  in  Figure  8 for 
the  1 MT  surface  burst  on  hard  rock.  This  parameter,  which  is  not 
provided  by  the  Post  and  McGetchin  et  al . models  (References  1 and  2, 
respectively),  exhibits  a stronger  range  dependence  than  the  depth 
distribution.  The  values  below  unity  may  be  interpreted  as  impact 
probabilities  (for  example,  there  is  a 10%  probability  of  one  impact 
per  m at  a range  of  1 km). 

The  primary  reason  for  the  high  values  of  the  areal  number 
density  in  the  crater  lip  region  (Figure  8)  is  the  preponderance  of 
small  fragments  at  the  short  ranges.  The  fragment  size  class  (defined 
by  the  minimum  and  maximum  size  fragments  reaching  a given  range) 
decreases  with  increasing  range  as  indicated  in  Figure  9.  This 
reduction  is  accompanied  by  a rapid  reduction  in  the  areal  number 
density  (Figure  8).  The  corresponding  minimum  and  maximum  values  of 
the  impact  time,  velocity  and  angle  are  shown  in  Figures  10,  11  and  12 
respectively. 
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Figure  11.  Ejecta  Impact  Velocity  for  i MT  Surface  Burst 
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Ejecta  cloud  outlines  for  the  1 MT  surface  burst  on  hard  rock 
are  shown  in  Figure  13  for  times  from  5 through  50  seconds  after  burst. 

The  ejecta  is  contained  in  a growing  fan-shaped  region  (in  two  dimensions) 
at  5 seconds  after  burst.  This  region  grows  but  maintains  a similar 
shape  to  about  20  seconds  after  burst.  After  that  time,  the  stem  begins 
to  form  along  the  vertical  axis  as  the  smaller  fragments  are  entrained 
by  the  rising  fireball.  The  larger  fragments  fall  toward  the  surface 
and  away  from  the  burst  point.  The  ejecta  cloud  outline  illustrates 
this  behavior  progressively  from  30  to  50  seconds  after  burst.  The 
entire  outer  region  of  ejecta  fragments  disappears  by  55  seconds  after 
burst,  leaving  the  stem  region  in  the  center. 

Density  contours  (number  of  fragments  per  m ) may  be  constructed 
from  the  ejecta  cloud  description  provided  by  the  current  model.  Examples 
of  these  contours  are  shown  in  Reference  11. 
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Section  6 


CONCLUDING  REMARKS 

A dynamic  model  for  the  high  explosive  and  nuclear  ejecta 
environments  has  been  presented,  together  with  selected  results  for 
a 1 MT  nuclear  surface  burst  on  hard  rock.  The  model  provides  both 
the  ground  deposition  and  airborne  ejecta  cloud  environments.  Predictions 
of  the  ejecta  depth  variation  with  range  compare  favorably  with  empirical 
correlations  of  other  investigators. 

The  ejecta  source  model  combines  current  theoretical  and 
empirical  results  in  a manner  that  permits  calculations  for  general 
cases  for  which  no  data  are  available  as  well  as  for  specific  high 
explosive  and  nuclear  events.  The  effects  of  explosive  type,  yield 
and  geology  are  taken  into  account. 

Two  specific  areas  where  further  research  would  lead  to 
significant  advancements  in  the  ejecta  model  have  been  identified. 

The  first  concerns  the  choice  of  source  parameters  for  bursts  in  incohesive 
soils.  The  use  of  the  in  situ  site  distribution,  including  the 
maximum  fragment  size,  will  give  substantially  different  results  than 
if  the  correlations  for  cohesive  soils  are  employed.  This  question 
requires  further  analysis  of  test  data  obtained  from  shots  over  dry 
incohesive  soils  for  its  satisfactory  resolution.  The  second  area  of 
concern  is  that  of  ejecta  rays.  The  scatter  of  nearly  four  orders  of 
magnitude  exhibited  by  models  based  on  correlations  of  measured  ejecta 
blanket  thicknesses  includes  azimuthal  variations.  The  current  ejecta 
model,  with  an  estimated  uncertainty  band  of  less  than  an  order  of 
magnitude, does  not  address  the  ray  problem.  Further  research  is  required 
to  determine  the  effects  of  rays  on  this  uncertainty  level. 
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NOMENCLATURE 


Symbol s 


a 

C 

Cr 


m 

M 

R 

Re 

SHOB 

T 

V 

W 

Z 

e 

y 

p 


Subscripts 

a 

al 

e 

m,  max 

min 

z 


Fragment  diameter;  term  in  Equation  (6) 
Coefficient  in  Equation  (12) 

Drag  coefficient 
Geology  correction  factor 
Source  correction  factor 
Mass  of  fragment 
Mass;  Mach  number 
Range;  radius 
Reynolds  number 
Scaled  height  of  burst 
Time 

Velocity 

Yield 

Altitude 

Angle 

Dynamic  viscosity 
Material  density 


Apparent  crater 

Apparent  crater  lip 

Ejecta 

Maximum 

Minimum 

Vertical 
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